Ultrasonic attenuation in clean d-wave superconductors. 
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We consider the attenuation of longitudinal ultrasonic waves in a clean two-dimensional el- 
wave superconductor. We show that the attenuation coefficient is linear in temperature at low 
temperatures for all in-plane directions of the propagation of the ultrasound, and that the coefficient 
of the linear term can be used to determine the parameters crucial for the low temperature transport 
in these compounds. 



Much of the novel physics associated with the super- 
conducting state of the high-T" c cuprates is due to non- 
trivial angular dependence of the energy gap, Ak, at 
the Fermi surface. While in conventional, s-wave, su- 
perconductors the energy gap is finite for all quasipar- 
ticle excitations, it is believed that in many cuprates it 
has the angular structure of a d x 2_ y 2 state, with lines of 
nodes, which leads to a gapless excitation spectrum along 
certain directions in momentum space. Consequently 
the temperature dependence of the thermodynamic and 
transport coefficients in the high-T c materials is quali- 
tatively different from that of their s-wave counterparts: 
at low temperatures the behavior of the coefficients mea- 
suring an average of a particular quantity over the Fermi 
surface is governed by the low-energy excitations near 
the nodes of the gap at the Fermi surface, and the re- 
sulting temperature dependence is given by power laws 
rather than by exponentially decaying functions with an 
activation barrierlj. While the power laws are entirely de- 
termined by the dimension of the manifold where the gap 
vanishes, in this case by the existence of lines of nodes, 
the coefficients can be used to determine parameters of 
the superconducting materials and test the agreement 
with specific models. 

The ratio of the Fermi velocity at the nodes, Vf, and 
the velocity vi associated with the growth of the super- 
conducting gap at the Fermi surface A(p/) = V2(p/ — 
PjT 1 ^) near the node is a particularly important parame- 
ter: both the universal limit of transport coefficientsi and 
the temperature dependence of the penetration depthel 
depend solely on v^/vf. The numerical value of this r, 
tio has not been clearly determined yet; Lee and We 
obtain Vf/v2 =ri6.8 from the analysis of the penetra- 
tion depth datacL while the measurements of the ther- 
mal conductivity!! yield Vf/v2 — 13.6. It is therefore 
important to have additional ways of experimentally de- 
tcrminining this parameter. 

Here we analyze the electronic attenuation of the lon- 
gitudinal ultrasound waves in a d-wave superconductor 
and argue that it can be used to measure the ratio Vfjv-i- 
We consider the clean case, ql 3> 1, where q is the 
wave vector of the sound wave, and I is the electron 
mean free path, this implies that for a linear frequency 
150MHz, taking the sound velocity v s = 4 x 10 5 cm/secl 
we need I > Afim. In previously studied optimally doped 



YBCO samples the mean free path was determined to 
be I ~ 0.6/im at T ~ 2OK0, which is closer to the hy- 
drodynamic limit, ql <C 1, where the longitudinal and 
transverse ultrasonic attenuation have been analyzedHB. 
However, the ultra .high purity crystals recently grown 
in BaZrOa cruciblesE3 are about an order of magnitude 
cleaner, and experiments demonstrate that at low tem- 
peratures I ~ 4 — b/irrEl, indicating that the limit ql > 1 
is likely to be achieved. Experimentally the clean limit is 
manifested in the attenuation rate which is linear, rather 
than quadratic, in the frequency of the ultrasound waves, 
and by an angular dependence of the longitudinal atten- 
uation rate, described below, which differs from that pre- 
dicted in the dirty limiC. 

Let us first consider the problem qualitatively. In the 
clean limit a s is determined by a scattering rate of a 
phonon by the quasiparticles. Apart from the matrix el- 
ement of the electron-phonon interaction this scattering 
rate depends solely on the phase space available for the 
phonon to decay. As the energy and the wave vector of 
the sound wave are small on the scale of electronic exci- 
tations, the relevant scattering processes are those with 
negligible energy and momentum transfer, and are given 
by the condition v ff ■ q = 0, where v g is the group velocity 
of the quasiparticles. In a normal metal the scattering 
is restricted to a belt on the Fermi surface where the 
Fermi velocity, v/, is perpendicular to q. The opening 
of an s-wave gap preserves the direction of the group ve- 
locity of the quasiparticles but thermally suppresses the 
occupancy of the states near the Fermi surface resulting 
in the exponential suppression of the scattering^. For 
an anisotropic s-wave superconductor the attenuation at 
temperatures lower than the minimal gap still decays ex- 
ponentially, ribeit with an exponent which depends on 
the directior£3. 

Recently this appmach has been transferred to a d- 
wave superconductoroEj. Assuming that the quasipar- 
ticles contributing to the attenuation are still located at 
the points of the Fermi surface where vt ■ q = 0, the au- 
thors of Refs. [l4|,|l5] noticed that the scattering processes 
sample the local gap in the direction perpendicular to q, 
and predicted the fourfold oscillations of a s as a function 
of the angle of propagation in the plane, 9. In particu- 
lar, since the gap vanishes along the nodal direction, the 
attenuation in that direction was predicted to be tem- 
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perature independent and equal to that at T c , while the 
attenuation in other directions was determined to decay 
exponentially with the activation energy given by the lo- 
cal gap A(0 ±7r/2). 

However for an anisotropic order parameter not only 
the magnitude but also the direction of the group veloc- 
ity changes with the opening of the energy gap. This 
change has important consequences for superconductors 
with gap nodes, especially for compounds with a rela- 
tively large value of the superconducting gap amplitude, 
A , as a fraction of the Fermi energy ej, and at tempera- 
tures below Ag/e/. This was first noted in the context of 
the heaw fermion superconductors by Coppersmith and 
Klemmta, who predicted a maximum in the ultrasonic 
attenuation for certain directions of the propagation of 
the ultrasound, and a power law behavior for low tem- 
peratures. Heavy fermion materials are not in the clean 
limit, and T c ~ IK made the desired regime impossible 
to achieve experimentally. In the cuprates, on the other 
hand, a high transition temperature and a large ratio of 
A /T c ~2.5-4 bring the gap amplitude to about 10-30% 
of the Fermi energy, so that the electronic contribution to 
the attenuation rate is measured at T <C Ag/e/, where 
the theoretical analysis for a d-wave superconductor is 
still lacking. 

We first proceed with the qualitative analysis. Be- 
low T c the excitation spectrum has the form E(k) — 
\ZCt + A k , where £k is the quasiparticle energy in the 
normal state with respect to e/, so that the direction 
of the group velocity, given by a tangential line to a 
E(k) = const contour, does not coincide with the direc- 
tion of the Fermi velocity at the same point. For T <C T c 
only the contours with E <C Ao, near the gap nodes, are 
important for scattering. Then, as shown in Fig. la, 
the contour of E = const is approximately an ellipse 
= ±^J E 2 — A' 2 ((j) — (j> n ) 2 , where </> parameterizes 
the Fermi surface, and A' is the angular derivative of the 
gap at the node, and for any wave vector q there are two 
quasiparticle scattering processes which contribute to the 
attenuation. Since a s is proportional to the phase space 
volume available for scattering an immediate conclusion 
is that the attenuation is linear in T at low temperatures 
in all the directions except normal to the nodes. 

A second observation concerns the role of the relatively 
large ratio Ao/ey. This can be accounted for by intro- 
ducing the curvature of the Fermi surface, as shown in 
Fig. lb. Then the group velocity of the quasiparticles 
with energies E <C A or E ^> A does not change sig- 
nificantly. For energies E ~ Ao, however, the competi- 
tion between the curvature of the contour E = const and 
that of the underlying Fermi surface opens up additional 
phase space for scattering of the ultrasound propagating 
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FIG. 1. a) A typical contour of constant energy near a 
node of the order parameter and the scattering processes con- 
tributing to the attenuation for an arbitrary q; b) Contours 
of constant energy E = aAo for a — 0.2; 0.8; 1.0; 1.4 with 
Ao/e/ = 0.1, and the scattering processes for q nearly nor- 
mal to the node. Arrows are parallel. 

within a narrow window of angles nearly perpendicular to 
the node; as seen in Fig. lb three rather than two scatter- 
ing processes become possible. For a d-wave gap the cur- 
vatures become comparable for E ~ A 2 /2ef ~ 2Ag/e/, 
which implies that the additional scattering is most ef- 
fective for T ~ 2A§/e/, a relatively high temperature, 
resulting in an attenuation coefficient larger than that 
at T c . For T < T c the additional scattering is effective, 
while at lower temperatures the occupancy of these states 
is thermally suppressed and they do not contribute to 
the attenuation. We therefore expect a maximum in a s , 
below the superconducting transition followed by a tran- 
sition to linear, in T, behavior at lower temperature, T*. 
The position of the maximum, as well as crossover tem- 
perature T* scale with Ag/e/, while the angular window 
within which it exists scales as Ag/e^. 

This analysis implies that the main conclusions of Ref. 
|i~o] are qualitatively applicable to d-wave superconduc- 
tors. We now investigate the issue in more detail focusing 
on the low-temperature regime. Within a standacd-BCS- 
like theory the ultrasonic attenuation is given byMEZl 



a s (T,q) = A 



n 

4T 



dk- 



coslr (£ k /2T) 



CkCk+q - A k A k+q 
-Ek-Ek+q 



(1) 



2 



0.4 




FIG. 2. Angular dependence of the low temperature slope 
of the attenuation for v^/tif = 0.2. 

where the constant A depends on the scattering matrix 



element and the frequency of the ultrasound, fl, and the 
integration is over the Brillouin zone. Expanding in q 
and taking into account that the angular average of the 
gap derivative vanishes, we obtain 

^-^h^kimi^-^ <2) 

where the group velocity v 9 is given by 

Vff (k) = V k £ k =-^(CkV / +A u ^). (3) 

We now consider a model two-dimensional d-wave super- 
conductor with the wave vector q in the plane at an angle 
9 from the x-axis. For q not exactly in the nodal or anti- 
nodal direction v/ • q = does not satisfy the condition 
imposed by the delta function; we exclude these two spe- 
cific cases from the analysis below, although we allow 9 
to be arbitrarily close to these directions, so that the re- 
striction has no bearing on the generality of the results. 
In the low-temperature regime T <C T c 



a(T c ] 



In 2 



v 2 T 



sin 2 (6»-7r/4) 



cos 2 (6»-tt/4) 



[cos 2 ((9 - tt/4) + (v 2 /v f ) 2 sm 2 {9 - tt/4)] 3/2 [sin 2 ((9 - tt/4) + {v 2 /v f Y cos 2 ((9 - tt/4)] 



3/2 



(4) 



Even though there is a normalization constant a(T c ) 
here, the ratio of the attenuation coefficients at different 
angles is a function of v 2 /vj only and therefore provides 
a direct measurement of this important parameter. No- 
tice also that the slope vanishes as v 2 /vf — > 0. Eq.fjJ) 
represents the main result of this Communication, and 
in Fig. 2 we show the slope as a function of the direction 
in the plane. 

To obtain an explicit temperature dependence and 
compare our results with those of Ref. 16 in more de- 



tail we now consider a superconductor with a cylindrical 
Fermi surface. Then the attenuation is given by 



sy > 2Tv f qJ 



i9-cj>)\ E((f>) 



cosh 



-2(E(4>)\ 
\ 2T )' 

(5) 



where Co(0) 



Agsin4<?itan(6' - (f>)/2e f , and E{4>) = 
VsoW + A 2 (0). Notice that the contribution to the 
integral of the regions cj) « 9 ± 7r/2, where vy ■ q = 
and cos(# — 0) tends to zero, is suppressed exponentially 
by the thermal function as tan(# — <p), and consequently 
E((f)), becomes large. 

In Fig. 3 we present the results for the ultrasonic at- 
tenuation normalized by its value at the transition tem- 
perature for near-nodal (as the nodes are at ra7r/4, the 



direction normal to a node coincides with the nodal di- 
rection) and anti-nodal orientation of the wave vector q 
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FIG. 3. Temperature dependence of the attenuation 
coefficient in the near-nodal (6 = 0.78, upper curves) 
and anti-nodal (8 = 1.57, lower curves) directions for ra- 
tios of Ao/ef of 0.25 (solid line), 0.1 (dashed line), 0.025 
(long-dashed line), and 0.0025 (dot-dashed line). 
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for different values of the ratio of Ao/e/ , which are very 
similar in general features to those obtained numerically 
in Ref. [ll| In the numerical work we have used the BCS 
ratio Ao = 2.14T C , a larger coefficient would push both 
the maximum and the onset of linear behavior towards 
higher temperatures. For the near nodal direction there 
is a clearly defined maximum of the attenuation at high 
temperatures, the position of the maximum scales with 
Aq/c/ as predicted above. The attenuation is linear at 
low temperatures. For the anti-nodal direction the decay 
of a s (T) is qualitatively close to the exponential behav- 
ior, although the dependence on the ratio Ao/e/ is clearly 
seen, and the shoulder on the curve for the largest value 
of the ratio indicates that the additional phase space for 
scattering has become available even for q near the anti- 
node. For both angles the results for the smallest ratio 
of Ao/e/ = 0.0025 are indistinguishable on the,seale of 
this graph from the generalized BCS predictiontL3~t3 that 
a s (T,0)/a n = /(A(0 + tt/2)), where / is the Fermi func- 
tion. We also obtain that for a fixed value of the ratio 
Ao/e/, the window within which the maximum can be 
observed is narrowly centered around the nodal direc- 
tion, as expected. Finally, Fig. 4 demonstrates the lin- 
ear, in temperature, behavior predicted here for T <C T c , 
which is in a sharp contrast to the exponential decay of 
the Fermi function, and agrees remarkably well with the 
result given in Eq.Q). 

The maximum of a s is at high temperature where the 
electronic contribution to the attenuation is difficult to 
measure. While our results are qualitatively the same 
for a different geometry of the Fermi surface, such as 
a tight-binding, the position of the maximum may be 
shifted; we also note that for a tight-binding Fermi sur- 
face close to half- filling the peak does not exist: large 
fiat regions of the Fermi surface away from the nodes 
contribute to the attenuation in the near-nodal direction 
above T c , and in the superconducting state the loss of 
phase space in the gapped regions cannot be compen- 
sated for by the increase in scattering near the nodes. 
Moderate impurity scattering is expected to "average" 
the attenuation over a small range of angles. As seen 
from Fig. 2, this can change the slope near the node, but 
has little effect on the slope away form the node, and 
this result should be robust with respect to scattering by 
dilute impurity centers. Therefore measurement of the 
low temperature slope provides a direct measurement of 
the parameter vijvf, which determines the low temper- 
ature behavior of the cuprates. We also note that we 
expect the general considerations of this work, including 
the breaking of particle-hole symmetry illustrated in Fig. 
Ib, to be important for other transport coefficients, such 
as Hall effect, although a different and detailed analysis 
is needed there. 
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FIG. 4. Low temperature dependence of the attenuation 
coefficient in the near-nodal (9 = 0.78) and anti-nodal 
(6* = 1.57) directions for Ao/e/ = 0.1. Dashed line: analytic 
result from Eq. (0) . 
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